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Abstract 

We study the nonrelativistic quantum Coulomb hamiltonian (i.e., 
inverse of distance potential) in R™, n = 1, 2, 3. We characterize their 
self-adjoint extensions and, in the unidimensional case, present a dis- 
cussion of controversies in the literature, particularly the question of 
the permeability of the origin. Potentials given by fundamental solu- 
tions of Laplace equation are also briefly considered. 

1 Introduction 

In principle the unidimensional (ID) hydrogen atom is a simplification of 
the three-dimensional (3D) model which has been invoked in theoretical and 
numerical studies [H El [3] ; note that Cole and Cohen [3] and Wong et al. 
[5] have reported some experimental evidence for the ID hydrogen atom. 
In a particular situation the ID eigenvalues coincide with the well-known 
eigenvalues of the 3D hydrogen model, as discussed ahead. 

Apparently the ID hydrogen atom was first considered in 1928 by Vrkl- 
jan [6]. However, it was a work of Loudon [7] published in 1959, whose 
potential model is 

Vo(x) = ~\~~\ ' K > °> 
\x\ 

that increased attention to the subject which has become interesting and 
quite controversial. We refer to Vc as the Coulomb potential. 

Loudon stated that the ID hydrogen atom was twofold degenerate, hav- 
ing even and odd eigenfunctions for each eigenvalue, except for the (even) 
ground state having infinite binding energy. Typically ID systems have 
no degenerate eigenvalues, and Loudon justified the double degeneracy as a 
consequence of the singular atomic potential. Andrews [8] questioned the ex- 
istence of a ground state with infinite binding energy. Ten years later Haines 
and Roberts [91 revised Loudon's work and obtained that their even wave 
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functions, with continuous eigenvalues, were complementary to odd func- 
tions, but such results were criticized by Andrews [TO], who did not accept 
the continuous eigenvalues. Gomes and Zimerman [Tl\ argued that the even 
states with finite energy should be excluded. Spector and Lee [12] presented 
a relativistic treatment that removed the problem of infinite binding energy 
of the ground state. Several other works P3J HH H51 HQ HZl OBI HSl 1201 121] 
(see also references therein) have discussed this apparently simple problem. 

In this work we advocate that the roots of such controversies is a lack of 
sufficient mathematical care in some papers: in ID the Coulomb singularity 
is so severe that it is not a trivial problem to assign boundary conditions 
at the origin. The main question is how to properly define the self-adjoint 
realization(s) of 

H = -^A + V c (x), domH = C$°(R\{0}). 

The domain choice is because functions ip in C^°(M\{0}) are kept far enough 
from the origin (i.e., zero does not belong to their support), and so Hip is 
well defined. H is hermitian but not self-adjoint, and it turns out that it 
has deficiency index n + = 2 = n_ and so an infinite family of self-adjoint 
extensions (see Section H} . Although such extensions appear in [TH] , details 
of how they were obtained are missing; in Section 0] we find such extensions 
by another approach, that is, we use a (modified) boundary form as discussed 
in [22] . These extensions are the candidates for the energy operator of the 
ID hydrogen atom. 

With these extensions at hand, we discuss the question of permeability 
of the origin, that is, whether in ID the Coulomb singularity acts as barrier 
that allows the electron to pass through it or not. This is one of the impor- 
tant questions considered in the literature. It is found that the permeability 
depends on the self-adjoint extension and we present explicit examples of 
both behaviors. 

In Sections [5] and we present the self-adjoint extensions of H in 3 and 
2 dimensions, respectively. It is well-known that in 3D the operator H, with 
the same action as H but domain C<^°(R 3 ), is essentially self-adjoint, that 
is, it has just one self-adjoint extension, whose domain is the Sobolev space 
W 2 (IR 3 ) (this is known as Kato-Rellich Theorem; see ahead). However, if 
the origin is removed and the initial domain C^°(M 3 \ {0}) is considered, 
then also in 3D there are infinitely many self-adjoint extensions. Note that 
in both ID and 2D the origin must be removed in order to get well-defined 
initial operators H. 

It is sometimes assumed that the right potential describing the coulom- 
bian interaction is given by the fundamental solutions of Laplace equation, 
that is, 

V\{x) = k\x\, V2(x) = rein \x\, Vs(x) = — -, — r, 
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in ID, 2D and 3D, respectively. For example, in the statistical mechanics 
of the Coulomb gas in 2D the potential V2 is often considered, instead of V3, 
and the so-called Kosterlitz-Thouless transition is obtained. So in Section [5] 
we consider the Schrodinger operator with such potentials and argue that 
they are always essentially self-adjoint (in suitable domains), independently 
of dimension. Finally, the conclusions are reported in Section [H 

A notational detail: the dot in H means that the origin has ben removed 
from the domain of the initial hermitian operator; e.g., in ID domH is 
Cg°(R\ {0}) and so on. 

2 Self-adjoint extensions: 3D 

The initial hermitian operator modelling the nonrelativistic quantum 3D 
hydrogen atom is (write r = \x\,6,ip for the spherical coordinates) 

h 2 

H = A + V c (x), domff = a 00 (M 3 )cL 2 (M 3 ), 

2m 

which is well defined since for if) S domH 

\\Vcii\\ 2 = k 2 [ ^ X }} 2 dx = k 2 f°° dr f dd sm6 d<p\ij;(r,9,<p)\ 2 < 00. 
Jr 3 Fl Jo Jo Jo 

The Kato-Rellich Theorem |23t [24"] applies to this case, since Vc £ 
L 2 (M 3 ) + L°°(M 3 ), and H has just one self-adjoint extension whose domain 
is the Sobolev space TC 2 (M?). Thus, the Schrodinger operator is well es- 
tablished in this case, so the quantum dynamics, and this is the standard 
operator discussed in textbooks on quantum mechanics (usually with less 
mathematical details). 

It is worth mentioning that in M. n , n > 4, Kato-Rellich Theorem implies 
unique self-adjointness for potentials V G L p (R n ) + L°°(lR n ) with p > n/2; 
so in dimensions n > 4 the Schrodinger operators H with potential Vc and 
domain Cg°(M n ) are always essentially self-adjoint. 

However, in ID and 2D the condition || Vc"0l| 2 < °o for all ip £ C^°(R 3 ) 
requires ip(0) = for x in a neighbourhood of the origin, that is, ip £ 
C^°(M n \ {0}), n = 1, 2. The self-adjoint extensions in such cases will be dis- 
cussed in other sections; for question of comparison with other dimensions, 
now we consider what happens if the origin is removed also in M 3 , that is, 
if the initial hermitian operator is 

H = -lf-A + Vc(x), domH = C °°(1R 3 \ {0}). 
2m 

Write £ = (9, ip) for the angular variables in the unit sphere S 2 and 
dt; = sin 9 d9 dip. Let V denote the set of linear combinations of products 
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f(r)w(0 with / G C£°(0,oo) C L 2 ((0,oo), r 2 dr) and w G L 2 (5 2 ,^). Due 
to the decomposition in spherical coordinates 

L 2 (M 3 ) = L 2 ((0, oo), r 2 dr) L 2 (S 2 ,d£), 

. For functions 4>(r,£) = f(r)w(£) G V the operator 



T> is a dense set in L 2 
i7 takes the form 



Hf(r)w(0 



£-1% + -%) f(r) - -f(r) 
2m \ r } r 



w{£) + 



2m r 2 



where B is the Laplace-Beltrami operator |25] 

1 



sin l 



dg (sin 6dew) + 



w 



sml 



acting in L(S 2 ,d£). B with domain C X (S 2 ) is essentially self-adjoint and 
its eigenfunctions are the spherical harmonics 5^ m (£)) which constitute an 
orthonormal basis of L 2 (5 2 ,d£); recall that 

(BY l>m )(^) = 1(1 + W, m (0, leN,-l<m<l. 

Denote by J\ the subspace spanned by {Y^ m : —I < m < I}, that is, the 
subspace corresponding to the eigenvalue 1(1+1) and := L 2 ((0, oo), r 2 dr)<& 
Ji; thus 

oo 

L 2 (M 3 ) = L 2 ((0,oo),r 2 dr)®L 2 (5 2 ,^) = 0Li- 

1=0 

If // is the identity operator on Si, the restriction of H to T>i = DnL; is 



given by H 



Hi <S> Ii, with 



H, 



2m V dr 2 



2 d 
r dr 



1(1 + 1) 



K 

r 



and our task is reduced to finding the self-adjoint extensions of Hi with 
domain Cq°(0, oo). It is convenient to introduce the unitary transformation 
U : L 2 ((0,oo), r 2 dr) -> L 2 (0,oo), (U(p)(r) = r</>(r), which maps Co°(0, oo) 
to itself and 

h 2 ( d 2 



hi 



UHiU- 1 



2m V dr 2 



1(1 + 1) 



K 

r 



with domain dom Hi = Cq°(0,oo). Standard arguments gives that the ad- 
joint h\ has the same action as h\ but with domain 

dom h\ = {<f) G L 2 (0,oo) : <p,<p' £ AC(0, oo), h*<f> G L 2 (0,oo)} . 

If f2 is an open subset of R, AC(f2) indicates the set of absolutely continuous 
functions in every bounded and closed subinterval of f2. 

By adapting the analysis of the free hamiltonian in M 3 with the origin 
removed, which was performed in [26], one proves the following result: 
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Theorem 1. ho has deficiency indices equal to 1, while hi, I ^ 0, is essen- 
tially self-adjoint. 

Hence, for I ^ the unique self-adjoint extension of hi is h*, while ho 
has infinitely many self-adjoint extensions. In order to find such extensions 
in case I = 0, we will make use of the following lemma \27\ [28l 129] . whose 
proof we adapt and reproduce. 

Lemma 1. If (j) £ dom/iQ, then the lateral limits 4>(0 + ) := lim r _^o+ 4>(r) and 

0(O + ) := lim ( 0'(r) H — 0(r) ln(ftr) 
r-*o+ V "> 

exist (and are finite). 

Proof. For £ dom h$ one has 

, * . h 2 d 2 (j) K . _ 2 . 

and one can write = 0i+02 with ^0" = u, 0i(O + ) = and 2"ir02+ K / r = 
0. Since <f>j G W 2 (e,oo), j = 1,2, for all e > 0, and u G L 2 , it follows that 
these functions are of class C 1 (0,oo). 

Consider an interval [r, c], < r < c < oo. Since 

0'i ( r ) - 0'i ( c ) = - h rJ r u ( s ) ds > 

(f>i(r) has a lateral limit 

1m f c 

/ 1 (o + ) = 0'( C ) + ^ y o u(«)d S . 

On integrating successively twice over the interval [r, c] one gets 
2 (c) -0 2 (r) = p- / — <k, 



and then 



02 (r) = </>2(c) - (c - r)0' 2 (c) - / dv [ ds^ 



V 



S 



= 02(c) - (c - r)0 2 (c) ^- y ds(f)(s) ——, 

and since < (s — r)/s < 1, by dominate convergence the last integral 
converges to Jq 0(s) as r — » + . Therefore 02 (0 + ) exists and 

02(0+) = 2 (C) - C0 2 ( C ) " ^ / 0(S) 
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Now, 



\<h(r) ~ M0 + )\ < rW 2 {c)\ + ^ r Ma)\ ds + ^r I ds ■ 



h 2 J '^'^ n 2 

Taking into account that (/> is bounded, say \4>( r )\ < C, Vr, Cauchy-Schwarz 
in L 2 implies 

(s)\ds = / l\<j)(s)\ds < Cy/r, 



and so, for r small enough and fixing c = 1, 

f ds^- < C(c\lnc\ +r|lnr|) < Cy/r, 

Jr S 

for some constant C. Such inequalities imply 0(r) = </>(0 + ) + 0{y/r), and 
on substituting this into 

2rriK f 4>(s) 

— ds 



(recall that 4>i(0 + ) is finite) it is found that there is b so that, as r — > + , 

0'(r) = 0'(1) - ^^(0+) ln(Kr) + b + o(l); 

thus, the derivative cj>' has a a logarithmic divergence as r — > and the 
statement in the lemma also follows. □ 

For tp £ dom Hq integration by parts gives 
where 

h 2 / , \ 

T(ip,cp) := lim ( tp(r)(f)'(r) — ifi (r)(j){r) ) 

2m r —>o+ V / 

is called a boundary form for /io |22j . Note that although T(ip,<j>) is finite, the 
lateral limit <^'(0 + ) can diverge; however, by Lemma Q] it is readily verified 
that 

r(V, 0) = ~ (v(o + )^) - ^(o+Mo)) 

and now all lateral limits are finite. The self-adjoint extensions of ho are 
restrictions of /iq to suitable subspaces D so that dom Hq C D C dom /iq and 
T\ D = 0, that is, T(i/i, <j>) = for all <p,ip e D. 

We introduce the unidimensional vector spaces X = {pi(ip) := ip(0 + ) + 
#(0 + ) : ij) G dom/i^} and Y = {p 2 (V0 := ^(0 + ) - #(0+) : ^ E dom/i^} 
and note that 

-fc2-r(V>,0) = {pi(lp),Pl(<f>))x - {f>2(lp),P2((f>)}Y, 
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with inner products in X, Y, as indicated. Hence, the subspaces D for 
which r vanishes are related to maps that preserve inner products, that is, 
unitary maps from X to Y (see details in |22| ) , and since these vector spaces 
are unidimensional such maps are multiplication by the complex numbers 
e lS , < 9 < 2ir. Therefore, for each 9 a self-adjoint extension of ho is 
characterized by the functions tp G dom/iQ so that P2(ip) = e l8 pi(i/j), that is, 

(l-e ie )V(0 + ) = i(l + e ie )4>{0 + ). 

If 9 this condition reduces to 

1 _|_ piO 

^(0 + )=A^(0+), A = ^Gl, 

and writing A = oo in case 9 = 0, the desired self-adjoint extensions /ig are 
described by 

dom/i^ = {-0 6 domh* : V>(0+) = A^(0 + )| , A G KU {oo}, 

and liQil) = h^ip. The Dirichlet boundary condition corresponds to A = 0. 
With such results at hand, we have 

Theorem 2. The self-adjoint extensions of H in 3D are 

oo 

H x = (u- l h^U ® J ) (U^h^U ® Ii) , AelU {oo}. 

i=i 

This should be compared with the case without removing the origin, 
for which there is just one self-adjoint extension. The eigenvalue equation 
for the Dirichlet case A = can be exactly solved in terms of Whittaker 
functions, and the negative eigenvalues are 

K 2 m 1 

n = "W^' n = 1 ' 2 ' 3 '-"' 

each one with multiplicity n 2 . For A / the manipulations become more 
involved and numerical procedures must be employed to find roots of implicit 
functions, and so the eigenvalues. 



3 Self-adjoint extensions: 2D 

Although our main interest is in the ID and 3D cases, we say something 
about the Coulomb system in 2D. As already mentioned, the origin must be 
excluded from the domain of the Coulomb potential in 2D and the initial 
hermitian operator is 

H = ~A + V c (x), domH = C °°(M 2 \{0}). 
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To find its self-adjoint extensions, introduce polar coordinates (r, ip) so that 



L 



L 2 ((0, oo), rdr) ®L 2 (S l ,dip), 



(S 1 is the usual unit circle in R 2 ) and the set T> of linear combinations of 
the products f(r)g(<p), f G C§°(0, oo) and g G C^S 1 ), is dense in L 2 (M 2 ). 
Now 



H}(r)g(y) 



2m r 2 



where ,6 = <9 2 is the Laplace-Beltrami operator acting in L 2 (5 1 ,dy?). This 
operator with domain C yo (S 1 ) is essentially self-adjoint, its eigenvectors 
gi(tp) = e ll{p /y / 27r constitute an orthonormal basis of L 2 (5' 1 , dip) and 

(B<B)(¥>) = -l 2 gi{ip), I G Z. 



Let [gi] denote the subspace spanned by gi and Li = L 2 ((0, oo), rdr)®[gi]; 



thus 



and if is the identity operator on [<#], the restriction of H to T>\ = T> n L; 



is given by £T 



Hi® II, with 



2m 



9 2 + -s r 

r 



- 



K 

: 

r 



with domain Cg°(0, oo), and the question is to find the self-adjoint extensions 
of such restrictions. By using the unitary operator U : L 2 ((0, oo), rdr) — > 
L 2 (0, oo), (U(j))(r) = r 1 / 2 (/>(r), one has 

fc- tr * trI =-IK* + (i-' , )p)-F 

with dom/i; = Cq°(0, oo) (since this set is invariant under U). By standard 
results it follows that the adjoint h\ has the same action as hi but with 
domain 

dom/i? = {0G L 2 (0,oo) : 4>,<p' £ AC(0, oo), fy> G L 2 (0,oo)} . 

Theorem 3. TTie operators hi are essentially self-adjoint if, and only if, 
I 0, whereas ho has deficiency indices equal to one. 



Proof. Weyl's limit point-limit circle criterion will be used 

i(j), that is, 



Thus we 



consider the solutions of hj 



— < 
2m 



2m V 4 



1 



0. 



G dom ft| . 
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By writing p = ^jpr, q = -j?f and performing the change of variable y 
(— 4g) 1//2 r, one gets 

0" + 



- - r — + -- - 

4 J y 2 y 4 

with r = 4C7) 1 / 2 . This equation has two linearly independent solutions 
given by the Whittaker functions \30\ 131 j 

My) = M TjW (y) e My) = ^r,\i\(y), 

whose asymptotic behaviors as \y\ — > oo are 

01 (y) ~ ey' 2 {-y)- T e 2 (y) ~ eT^V- 

Since there is no c G M so that 02 G L 2 (c, oo), it follows that hi is limit 
point at oo. Note that the above asymptotic behaviors as \y\ — > oo do not 
depend on I; however, at the origin we need to separate the cases I = and 
I + 0. 

For I = [301 ED, 

M0 + ) = and 2 (O + ) = 0. 

In this case there is c > so that <j)\, <p% G L 2 (0, c) and /to is limit circle at 0. 
Therefore, ho is not essentially self-adjoint but has deficiency indices equal 
to 1. 

For I 7^ 0, (pi (0 + ) = while <fe(y) diverges as 

E 1 Tm k ^ k) nk -\l\~r+ l/2)(-y)-^+\ 

k=0 

for y — > + (here T denotes the well-known Gamma function). Therefore, 
/i; (i 7^ 0) is limit point at 0. By Weyl criterion, hi is essentially self-adjoint 
if I + 0. □ 

Since hi, I 7^ 0, is essentially self-adjoint, its unique self-adjoint extension 
is exactly hf . According to the proof of Theorem El the deficiency subspace 
K-(ho) is spanned by (f)-(r) = M. T fl{{— Aq) l l 2 r) and the deficiency subspace 
K + (ho) spanned by 4>+(r) = 0_(r). The von Neumann theory of self-adjoint 
extensions [321 [24"] characterizes them by unitary maps from, say, K- to 
K + , and since such subspaces are unidimensional these unitary maps are 
just multiplication by e %e , < 9 < 2tt. Thus there is a self-adjoint extension 
/ig of ho, for each 9, and if ho denotes the closure of ho, it is given by 

dom/iQ = j-0 + c - e i( V-) : V> G domT^, c€ c| , 
^oV = ^ l oV ; ) V> ^ dom/iQ. 

In summary: 
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Theorem 4. The self-adjoint extensions of H in 2D are 

oo 

H e = (u-^UQIq) (U-^-hfU ® Ii) , 0e[O,27r). 
4 Hydrogen atom: ID 

In the unidimensional case the initial hermitian operator is 

and the origin naturally decomposes the space 

L 2 (IR \ {0}) = L 2 (-oo, 0) © L 2 (0, oo), 

and also the domain of H into C^°(— oo, 0) and C^°(0, oo); let H+ and H- 
denote the restriction of H to these subspaces, respectively. Thus, we have 

H = H-®H + . 

From the physical point of view, an important question is about the behavior 
of the system at the origin, e.g., is it impermeable, so that the system 
actually decomposes into a right one and a left one, or is it permeable? In 
the latter possibility, what do happen with wavefunctions at the transition 
point (the origin)? As discussed ahead, there are plenty of possibilities, due 
to infinitely many self-adjoint extensions. 

4.1 Self-adjoint extensions 

The adjoint operator has the same action as H but with domain 

domij; = {(f>e L 2 (0,oo) :<f>,<f>' € AC(0,oo), H* + <j) G L 2 (0,oo)} , 

and an analogous expression for H*_ and its domain. 

To find the deficiency subspace K + (H + ) we look for solutions of 

that belong to domi/^. Write p = q = -^r and perform the change 
of variable y = (— 4g) 1//2 x, so that this equation takes the form 
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with t = p/i—^q) 1 ^ 2 , which has exactly two linearly independent solutions 

PH EH 

0i+ (y) = Wt,i/2(v), 02+(y) = M Til/2 {y)- 

These solutions have finite limites as y — > 0. The asymptotic behaviors for 
\y\ — > oo are 



1+ (2/)~y^/ 2 and fc+fc/) ~ (-y)~ T 



2//2 



Hence, in the original variable the deficiency subspace K + (H + ) is unidi- 
mensional and spanned by <fii+(x) = W T ,i/2{(~ ^q) l ^ 2 x). A similar analysis 
implies that K^{H + ) is also unidimensional and spanned by <p\ + (x). There- 
fore H + has both deficiency indices equal to 1. 

Similarly one finds that the deficiency subspaces K + (H_) and K_(H_) 
of H_ are spanned, respectively, by = VV T;1 / 2 ((— 4g) 1 / 2 |x|) and (f>i-(x). 

Therefore H- also has both deficiency indices equal to 1. 

Now, due to the above results, the adjoint operator H* has the same 
action as H but domain (write TL = L 2 (IR \ {0})) 

domH* = \ G H :<j>,<j>' G AC(M \ {0}), <f>" - -~4> G H 

y 2m \x\ 

and both deficiency subspaces K±(H) have dimension 2; K + (H) is spanned 
by 

^to-/ if x>0 and nbo(x\-[ if x > 

if x<0 and t2[x) ~ \ if x<0 ' 

and K-{H) by ?/>i and ip%. 

Therefore, H has both deficiency indices equal to 2 and it has infinitely 
many self-adjoint extensions. A boundary form will be used in order to get 
such extensions. The following lemma will be needed, whose proof is similar 
to the proof of Lemma Q] (note the difference of signs in the definitions of 
0(0+) and 4>(0~) below). 

Lemma 2. If (j) G domH*, then the lateral limits 0(O ± ) := lim^^Qi 4>{x) 
and 

0(0*) := lim <//(z) ± — =-0(x) ln(±Kx) 
exzsi and are finite. 



For ijj,<j) €L domH* one has, upon integrating by parts, 
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where 

-f?lW) 



lim ( ip(x)4>'(x) — ip'(x)(f)(x) ) + lim ( — ip(x)(j)'(x) + ip' (x)(j)(x' 
and, using Lemma [21 straightforward computation gives 



rOM) = (^(o + )<A(o+) - V(o + )^(o+) - ^(o-)0(o-) + ^(o-)^(o- 

and now each lateral limit is finite. 

Introduce two linear maps pi , P2 : dom iJ* — > C 2 : 

/ vJ(o+) + #(o+) \ j , n /^(o+)-#(o+) 

so that 

(pi(VO>Pi(0))c a - (P2(V0,P2O))c 2 = --^2"« r (V ; ,^), V^,0 e domH*. 

As in Section O the self-adjoint extensions of H are restrictions of -ff* to 
suitable subspaces D so that domif C D C domij* and r|/) = 0, that is, 
r(V>, 0) = for all cf>, ^ G D [22]. 

Vanishing of T on domains D is equivalent to the preservation of the 
inner products in C 2 , and so it corresponds to unitary 2x2 matrices U, 
and each of such matrices characterizes a self-adjoint extension of H, 
so that domi?^- is composed of ip £ domH* so that p2(ip) = Upi(ip); also 
Hjjip = H*ip for if) £ domHy. 

The condition P2(V0 = Upi(t/j) is then written 

and we have explicitly got the boundary conditions characterizing the de- 
sired self-adjoint extensions. 

In case (I — U) is invertible (similarly if (I+U) is invertible) it is possible 
to write the above boundary conditions in the form 



VK0-) J V ^(°~) 

with A = — ?7)~ 1 (J + U) being a 2 x 2 self-adjoint matrix. 

In [18] the authors have got this form for the self-adjoint extensions, and 
it was claimed that by allowing the entries of A taking infinity all self-adjoint 
extensions are found; we think it is a hard task to cover all possibilities above 
(i.e., via U) with this representation via self-adjoint matrices A. 
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Particular choices of the matrix U (I is the identity matrix) 



a)/, b)-/, c) (J J), d) (_» 

impose, respectively, the boundary conditions: a) V^O - ) = = tp(0 + ) 
(Dirichlet); b) tp(0 ) = = ^(0+) ("Neumann"); c) V(0") = ^(0+) and 
J(0~) = J(0+) (periodic); d) V(0") = -ip(0+) and ^((T) = -</>(0+) (an- 
tiperiodic). 

Recall that the general form of a 2 x 2 unitary matrix is 

U = e i8 (^ a b #e[0,2vr), a, b G C, |a| 2 + |6| 2 = 1. 

This form will be used ahead. 



4.2 Negative eigenvalues 

In this subsection we discuss the negative eigenvalues of some self-adjoint 
extensions Hjy. The main goal is to remark that the eigenvalues, their mul- 
tiplicities and parity of eigenfunctions depend on the boundary conditions. 
This becomes important since in the past some authors have assumed partic- 
ular hypotheses on the eigenfunctions (see references in the Introduction), 
but without specifying the self-adjoint extension they were working with; 
this was the main source of controversies in the studies of the unidimen- 
sional hydrogen atom. 

Our first analysis is for U = I, i.e., Dirichlet boundary condition, and 
denote by Hu the operator 

h 2 d 2 k 

H D = "^^2 - 0> domH D ={<j>e domH* : 0(0+) = = 0(0")}. 

As in [18] we consider the Green function of {Hj^ — E) -1 , denoted by G(x, y), 
that is 

(H D - £)-Vx) = 9(x) f G(x, y)u(y)dy + &(-x) f G(x, y)u(y)dy. 

JO Jx 

Here Q(x) = 1 if x > and vanishes if x < 0. 

If u G mg{H£> — E), we search for solutions <f> of 

(H D - E)4> = u, (1) 

and the discussion is for (— oo,0) and (0, oo) separately. 

For x G (0, oo), by the method of variation of parameters, one finds the 
solution 

<f)(x) = fa(x) / - , dy + <t> 2 (x) / — dy, 

Jo W x {(f>i,(t>2) Jo W x {<pl,4>2) 
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where <p\ and cp 2 are independent solutions of the homogeneous equation 
(H D - E)<j) = 0, that is, 

-|/-(; + *)*=°' < 2 > 

and W x (4>i,4>2) the corresponding Wronskian. 

Writing p = q = ^f, r = p/i-Aq) 1 / 2 and z = {-Aq) l / 2 x, the last 
equation takes the form 

whose independent solutions are 

M z ) = W Ttl/2 (z) and <p 2 {z) = M T>1/2 (z). 

Recall that W Tjl / 2 (z) ~ e~ z / 2 z T and -M Tjl / 2 (z) ~ e z / 2 (-z)~ T , as z — > 00. 
In the original variable 

^1,02) =4^> 

and the unique solution satisfying 0(O + ) = is 

= [ (Wr,i/2((-±q) l/2 x) M T , 1/2 {(-Aq) l l 2 y) - 

M Tjl/2 ((-4q)^ 2 x) W r ,i/ 2 ((-4g) 1/2 y)) u(y)di/. 



Similarly, for x € (— 00, 0), the unique solution satisfying 0(0 ) = is 

= [ ^V2((-49) 1/2 kl) -M.,i/ 2 ((-4g) 1/2 |y|) - 

M r ,i/ 2 ((-4g) 1 / 2 | ; r|) W r , 1/2 ((-4 (? ) 1 / 2 |y|)) ufojdy. 

Summing up, the Green function of the resolvent operator (Ho — E)^ 1 
is given by 

G(x,y) = e(xy)^l 
Q(\x\ - |y|)W Til/2 ((-4g) 1 / 2 |x|)^ ril/2 ((-4 (7 ) 1 /2| y |) _ {x „ y) 

The values E for which (Ho — E)^ 1 does not exist constitute the eigen- 
values of Ho, and they are obtained from the points for which the Gamma 
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function T(l — r) is not denned, that is, 1 — r is a negative integer num- 
ber. By recalling the expressions of p, q and r, the condition 1 — r = — n, 
n = 0, 1, 2, 3, ■ ■ -, gives 

K 2 m 1 

-E'n = o n = 1) 2, 3, • • • , 

2rr n z 

which coincide with the eigenvalues of the usual 3D hydrogen atom model. 

These eigenvalues are twofold degenerated and a basis {(ft n ,i, 4>n,2} of the 
subsequent eigenspace is 

Mx) = &((-l) k x) W r , 1/2 ((-4g) 1 / 2 |x|), k = 1, 2. 

The negative eigenvalues of other extensions Hjy are harder to get and 
numerical computation should be employed. In the following particular cases 
of interest are selected. Let ^S>(x) := ^(lnT(x)) and define 



2mn 



2 



:n(|-r) '-"I' I 'I' I - r)- I 



(-2 J Em) 1 /2 



ft 



with [r(l — t)] 1 and ±c*;(£?)[r(l — r)] 1 denoting, respectively, the values 
of the lateral limits lim^Qi W Tj i/2((— 4g , ) 1 / 2 |x|) and 

Jim (A Wr l/2(( _ 4(? )V2 W) ± ^ Wr l/2(( _ 4(z) i/2| x | )ln(±ra ^ . 

Given a unitary matrix [7, the candidates for eigenfunctions of Hjy must 
satisfy the corresponding boundary conditions. 

Example 1. Let's take 6 = |, a = 1 and 6 = so that 

it ( 1 

and the values of £7 for which u{E) = — 1 are the eigenvalues of Hy and 
with multiplicity two; the corresponding eigenspace is spanned by 

Mx) = @((-l) k x)W TA/2 ((-4q) 1 / 2 \x\), k = 1, 2. 

Example 2. Consider another case: = f , a = i and 6 = so that 



i 
-i 



and the values of E for which u{E) = are found to be nondegenerate 
eigenvalues of Hy, and for each eigenvalue the corresponding eigenspace is 
spanned by 

4>{x) = e(x)w Til/2 ((-4 9 ) 1 /2| a; |). 
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These two cases illustrate that the behavior of eigenfunctions are not 
related only to the parity of the potential Vc(x), since there are cases for 
which the eigenfunctions do not have a definite parity and cases with eigen- 
values simple as well as with multiplicity two. These different possibilities 
are directly related to the singularity of the potential and depend on the 
selected self-adjoint extension. 



4.3 Permeability of the origin 

Another question that has been discussed in the literature is about the 
permeability of the origin in the unidimensional hydrogen atom; see, for 
instance [H E31 EH]- Some authors consider the origin an 

impermeable barrier, while others assume it is permeable. Again the answer 
strongly depends on the self-adjoint extension considered (see also [55]). 
as illustrated ahead. Here the definition of permeability is through the 
probability current density; for simplicity we assume h = 1 and m = 1. 
Recall that the probability current density j(x) in ID is given by 

j(x) = | (4>{x)<j)'(x) - (f/(x)<j)(xyj , (j) £ domij^, 
and it satisfies the continuity equation 

§- t m,x)\ 2 + ^j(t,x) = o. 

Our previous results in Subsection 14.11 show that lim x ^ ± j(x) do exist (see 
ahead). 

For each <j) £ domHjy, on integrating by parts we get 

= (H^A) - (&H & <f>) = ihm[j(e)-j(-e)}. 

Hence, the function j(x) can be continuously defined at the origin j(0) via 
lateral limits. Physically this relation means that the current density is 
isotropic at the origin, in spite of the strong singularity there. A simple 
observation shows that 



j(x) = Im((j)'(x)cj>(x)), 4> G dom Hjy, 
where Im indicates imaginary part, so that 



j(0) = lim Im ((f> (x)(f)(x)) = lim Im (cj)' (x)(j)(x)). 

Since 4>'(0 + ) and 0'(O~) can be divergent, we use using Lemma [2] to 
obtain 

j(0) = lim Im((p(x)(j)(x)) = lim Im (<f>(x)(j)(x)). 

x^0+ x^0~ 
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Note that it is exactly this relation that guaranties that j(0) is well defined 
and finite. 

We are now in position of giving a rigorous definition of permeability: If 
j(0) = 0, £ dom Hjy, the electron is completely reflected when approach- 
ing the origin, and so we say the origin is not permeable (or is impermeable), 
so that the regions x < and x > are kept separated by the singularity. 
If j(0) / we say the origin is permeable. 

Next we study the current density related to Hy in three cases. 

Case 1. (I — U) is invertible. Since A = — L r ) _1 (/ + t r ) is a self-adjoint 
matrix, the boundary conditions of Hjj become 

0(0+) Wu z\(-4(o*)\ uveRzeC 

4>(o-) ) \z v ) y 0(o-) )> MeK,ZGL ' 

and u, z and v are functions of the entries of U. The boundary conditions 
become 

lim 4>(x) = —u lim Mx) + z lim Mx) (3) 

x^0+ x^0+ x->0~ 

lim 4>(x) = —~z lim (j)(x) + v lim 4>(x) (4) 
x^o- x^0+ x^0~ 



Multiply the first equation (before taking limits) by <p(x) to get 



lim <j)(x)(j)(x) = — u lim |0(x)| 2 + z lim (j)(—x)(f>(x), 

x^0+ x^0+ x^0+ 



and so 



j(0) = lim lm((f)(x)(f)(x)) = lim lm(z(f)(— x)(j)(x)). 

x->0+ x^0+ 



Therefore, if z = then j(0) = 0, V0 G domi^, and we have found 
a family of self-adjoint extensions for which the origin is not permeable. 
Example [1] above corresponds to the self-adjoint matrix 



-U)-\I + U) = 
and the origin is not permeable in this case. 



1 
1 



Case 2. (I + U) is invertible. The matrix A = i(I + U)~ l {I - U) is also 
self-adjoint and the boundary conditions of Hjj take the form 

u z\f 0(0+) \_( -0( o+ ) ^ uvGMzGC 

and u, z and v are functions of the entries of U. Table 1 shows the current 
density at the origin for various values of u, z and v. 
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Table 1: Current density for invertible (I + U). 



u, z, v 7^ 


j(0) = lm(-flim^ - faM-x)) 


z = 


j(0) = o 


z 7^ and u = 


j(0) = Im (-i lim a ._ >0 + ^(x)^(-a:) 


> 


z / and v = 


i(0) = Im (| linx^ - 0(x)0(-x)j 





Note that Dirichlet boundary condition (so U = I) is a particular case 
with z = (the matrix A = 0), and since this case can not be an extension 
of another self-adjoint extension of H, we conclude that the current density 
vanishes for all (f) 6 domii^ if, and only if, z = 0. In other words, if (/ + [/) 
is invertible, the origin is impermeable precisely if z = 0. As expected, 
Dirichlet boundary condition implies the origin is impermeable. 

Case 3. Both (I + U) and (I — U) are not invertible. This case amounts to 

det(/ + L>) = = det(/-f7), 
which turns out to be equivalent to the following matrix representation 

U = ( ~_y V ) , u£M,v£C, Id 2 + \v\ 2 = 1. 

\ V u J 

The current density always vanishes at the origin if, and only if, v = 0, that 
is, the matrix U equals 

/ -1 \ ( 1 

{ 1 J or (,0-1 

If v / 0, we have 

lim 6(x)(b(—x) — lim 6(x)(b(— x) + i lim \6{ 

1 + !1^0- 1 / l ' (1 + ll) x-O" ' l+Ux^0+* Yy 



Example 3. Consider the self-adjoint extension which corresponds to 
the unitary matrix 

y/2/2 -y/2/2 
y/2/2 V2/2 



X)! 2 



U = i 



The domain of Hjy constitutes of the ip £ dom H* so that 
^(0+) \ = (V2 -i\f -v(o+ 

^(o-) J \ i V2 J V ^(o-) 
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and such conditions imply 



j(0) = lim Im(— x)ip(x)). 

The values of E for which oj{E) = —y/2 are eigenvalues of H^y of multiplicity 
two, and the corresponding eigenfunctions are 

<j> k {x) = e((-l) fe z)W T , 1/2 ((-4g) 1/2 |x|), k = 1,2. 

By taking the linear combination ^(x) = (fii(x) +4>2(x), and the asymptotic 
behavior of such eigenfunctions near zero, discussed in Subsection 14.21 we 
obtain 

j(o) = -r(i-r)-Vo, 

that is, if the electron is in this eigenstate it is transmitted through the 
origin. 

We conclude that there are extensions for which the origin is permeable 
and for others it is impermeable. Andrews [5] defines j(x) = i[F'{x)F{x) — 
F(x)F'(x)] but computes only for eigenfunctions; it is clear that j van- 
ishes if the eigenvalue is nondegenerate, since the corresponding eigenfunc- 
tion can be taken real. Andrews mentioned the possibility of zero current 
in case of degenerated eigenvalues. In our Example [T] the eigenvalues have 
multiplicity two and the current density is zero, whereas for the operator in 
Example [3] the eigenvalues have multiplicity two and the origin is permeable; 
therefore both possibilities are allowed in case of multiple eigenvalues. 

We note that the analysis of Moshinsky [29| . although interesting, consid- 
ers the "eigenfunctions" W\(—z) and W\{z) that do no belong to L 2 (M\{0}). 

5 Potentials Via Laplace Equation 

This brief section deals with Schrodinger operators with potentials V along 
the fundamental solutions of Laplace equation 

AV = 0. 

As mentioned in the Introduction, in physics sometimes one assumes 
that in each dimension the potential describing the Coulomb interaction 
is the fundamental solutions of Laplace equation [35l [36] ; recall that these 
solutions are (take k > 0) 

Vi(x) = k\x\, V2(x) = rein \x\, V^(x) = —- — r , 

\x\ 

in ID, 2D and 3D, respectively. 
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The case of V3 in 3D is standard and was recalled in Section [2j we 
underline that the operator 

h 2 

H = -—A + V 3 (x), domtf = C CO (R 3 ), 
2m 

is essentially self-adjoint and its unique self-adjoint extension has the same 
action but domain 7Y 2 (IR 3 ) ; this extension has nonempty discrete and essen- 
tial spectra. 

The case of potential V2 was analyzed by Gesztesy and Pittner [37] and 
they state the following result: 

Theorem 5. The operator 

k 2 

H = -—A + Khx\x\, domH = C%°(R 2 ), 

is essentially self-adjoint, and its unique self-adjoint extension has empty 
essential spectrum. 

Now we consider the unidimensional Schrodinger operator 

h 2 d 2 



00 



H = — r + k\x\, domH = C 

2m dx z 

whose adjoint H* has the same action as H but domain 

domff* = {<f> G L 2 (R) : fatf e AC(R), H*(j) G L 2 (M)} . 
Since V% is a bounded from below and continuous potential, with 

lim V\{x) = 00, 

I as I — >oo 

the following theorem follows from general results [231 [23] . 

Theorem 6. The above unidimensional operator H is essentially self-adjoint, 
its unique self-adjoint extension H* is bounded from below and has empty 
essential spectrum. 

By solving the eigenvalue equation 

h 2 

0" + (k|z| -E)<f> = 

2m 

in terms of Airy functions we have found the eigenvalues E n are simple and 
with asymptotic behavior 



2m 



mn 3tt , , 



2/3 

n — > 00. 



Hence, for Schrodinger operators H with potentials along the fundamen- 
tal solutions of the Laplace equation, we have: 
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1. H is essentially self-adjoint in C^(1R 3 ) and its self-adjoint extension 
has both nonempty discrete and essential spectra. 

2. H is essentially self-adjoint in C§°(M. n ) and its self-adjoint extension 
has purely discrete spectrum for n = 1, 2. 

However, for Schrddinger operators whose potential is the Coulomb one, 
i.e., Vc(x), we have: 

1. The deficiency indices are equal to in C^°(1R 3 ). 

2. The deficiency indices are equal to 1, 1 and 2 in C£°(R 3 \{0}), C$°{R 2 \ 
{0}) and Cg°(R \ {0}), respectively. 

6 Conclusions 

Although the 3D usual model hamiltonian H with Coulomb potential Vc, 
domH = Cq°(M 3 ), is essentially self-adjoint, in W 1 the l/\x\ singularity 
imposes the initial domain must be C^°(IR ra \ {0}), n = 1, 2, and the corre- 
sponding operators H have deficiency indices equal to 2 an 1, respectively; 
hence with infinitely many self-adjoint extensions. For the sake of compar- 
ison, we have also considered the origin removed in R 3 , that is, H with 
domain C^°(R 3 \ {0}): the deficiency indices are equal to 1 in this case. In 
each case, all self-adjoint extensions have been found. 

In ID the question of permeability of the origin was analyzed and the 
answer depends strongly on the self-adjoint extension considered. Due to 
particular examples discussed, we conclude that the multiplicity two of the 
eigenvalues does not determine the permeability. 

We have paid particular attention to the ID case, since there are many 
papers in the literature about this model and occasionally with conflicting 
conclusions. We have found that these conflicting positions have been orig- 
inated from boundary conditions imposed mainly on "physical basis" that 
can fail for strong singularities, as is the case of Vc in one-dimension. We 
expect to have clarified the situation, and the next step could be present- 
ing arguments to select the extension (s) to be considered natural, with the 
consequent implications as, for instance, the permeability of the origin. 

Finally, we have found remarkable that, for potentials in R n , n = 1, 2, 3, 
given by fundamental solutions of Laplace equation, the corresponding initial 
hermitian operators with domain C^°(M n ) are always essentially self-adjoint. 
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